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A recently presented anisotropic generalization of the multicomponent super- 
symmetric t — J model in one dimension is investigated. This model of fermions 
with general spin-S is solved by Bethe ansatz for the ground state and the low-lying 
excitations. Due to the anisotropy of the interaction the model possesses 2S massive 
modes and one single gapless excitation. The physical properties indicate the exis- 
tence of Cooper-type multiplets of 25 + 1 fermions with finite binding energy. The 
critical behaviour is described by a c = 1 conformal field theory with continuously 
varying exponents depending on the particle density. There are two distinct regimes 
of the phase diagram with dominating density-density and multiplet-multiplet cor- 
relations, respectively. The effective mass of the charge carriers is calculated. In 
comparison to the limit of isotropic interactions the mass is strongly enhanced in 
general. 
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1 Introduction 



The t — J model is one of the most fundamental models of strongly correlated electron 
systems. It has been studied intensively in recent years since its connection with high 
Tc superconductivity was pointed out [ 1, 2]. It has been conjectured [ 3] that one 
and two-dimensional variants of this model have properties in common. Therefore it 
is important to have exact results available in one dimension for integrable coupling 
parameters. 

The t — J model describes electrons with nearest-neighbour hopping and spin- 
exchange interaction with the constraint that no site is allowed to be occupied by 
more than one particle. The one- dimensional model for spin-1/2 particles was found 
to be integrable at the supersymmetric point [ 4- 6] and the critical exponents 
for the long-distance asymptotics of its correlation functions were calculated [ 7]. 
These results were extended to a t — J model of fermions of arbitrary spin S with 
supersymmetric interaction [ 8- 10]. Recently, the integrability of an anisotropic 
generalization of the multicomponent supersymmetric t — J model has been shown [ 
11, 12]. The Hamiltonian reads 

H= - ^ P (c+ 9+1, + c++i,Cj,) P 

+ {4s^js'C'j+is'Cj+is - exp[sign(s' - s)r]]njsnj+is'} , (1) 

where Cjg annihilates a fermion at site j with spin component s. We assume s = 1, 
2, 25" -|- 1 for convenience, rijs denotes the number operator and V the projection 
operator onto the subspace of at most singly occupied sites. For 77 = the interaction 
is supersymmetric. In the general case the Hamiltonian is sp 1^(25" -|- 1, 1) invariant 
and ?7 > is a measure of the anisotropy. 

In this paper we calculate the critical exponents of the correlation functions for 
the integrable anisotropic t — J model of fermions with arbitrary spin 5* (|l|). For 
5* = 1/2 some of these results were presented in [ 13]. The paper is organized as 
follows. In section 2 the Bethe ansatz solution is presented and the ground state 
and excitations are discussed. In section 3 the density-density and superconducting 
correlation functions are investigated, and the transport masses are calculated. The 
paper ends with a discussion of the results. 

2 The Bethe Ansatz 

The general Bethe ansatz solution of the eigenvalue problem for Hamiltonian (|I]) 
has been obtained on the basis of the Perk-Schultz model [ 14- 16] in [11] from 
which we quote the relevant equations. The eigenstates are characterized by sets of 
rapidities A^^-* (j = 1, 2, A^) for the electrons which total number is N, and by spin 
rapidities A^'^ {I = 1, 2, 25*; a = 1, 2, Mi). The eigenvalue problem reduces to 
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the task of solving the following Bethe ansatz equations 
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where L is the number of lattice sites and we have set Mq = N and M2S+1 = 
0. (The number of particles with spin component m is conserved and equal to 
M2s+i-m — M2s+2-m-) The total energy and momentum of the system are given in 



terms of the electron rapidities A,- 
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where we have used the definition 
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We now discuss the patterns of distributions of rapidities satisfying (Q). Due to the 
attractive nature of the interaction the ground state is decribed by bound complexes 
of rapidities A^') (/ = 0, 1, 2^ - 1) which can be parametrized conveniently by 
the (real) rapidities A*-^'^-'. (In this paragraph subscripts are dropped). In order 
to accomodate for the groundstate we assume that the number of electrons is a 
multiple of 25* + 1, iV = M{2S + 1). Then we have M25 = M real rapidities A^^-^^ 
and for each such A*^^'^-' there are 25* + 1 — / complex spin rapidities A*^'^ given by 



A(2^) + ^? 
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-{23 - I), -{23 -1-2), {23-1- 2), (25 - /). (5) 



In particular we have Mi = M{23 + 1 — /) from which we conclude that the particles 
are equally distributed over all spin components, i.e. the state is a singlet. 

Using d^) the equations @ are straightforwardly reduced to one set of M equa- 
tions for the parameters Va = 2X^^^^ 



LQ [Va] — - — V 
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25 

e{v) := (6) 

1=1 

The numbers la are integer (half-integer) for even (odd) 2S{M — 1). In the ground- 
state these numbers are arranged symmetrically around zero. The energy and mo- 
mentum are given by 

M M 

E=Y.e{va), P=Y.p{va), (7) 

a=l a=l 

with the following abbreviations 

J sinh(25 + l)r/sinhr/ /oc , ia u \ 

e{v) := 2\ — . . (2^ + l)coshr7^ 

[ cosh(2D + Ijr/ — cosf J 

p(v) := TT - e (^v; 2^ ^^ 

In the thermodynamic limit L, M —>■ oo with finite M/L the parameters Va fill 
an interval [vq, 2it — Vq] uniformly with density cr{v). /^From (^) we obtain an integral 
equation of Fredholm type for the distribution function a{v) 

/ 2 9-1-1 \ n2Tr—vo 

2na{v) = Q' (v; —v)- 9'{v - w)aiw)dw, (9) 

V 2 J JvQ 

with Q'[v,ri) = sinh 2?7/(cosh 2// — cosf) and subsidiary condition 

where p is the particle density. The solution to these equations a{v) yields the ground 
state energy per site 

En r2-K-vo 

^ = / e{v)(T{v)dv. (11) 

A gapless mode of excitations is given by a redistribution of the Va parameters 
in (H) by keeping the string structure (^. These excitations are of particle-hole type 
with quasilinear dispersion. There are 25" massive modes of excitations corresponding 
to the breaking of some complexes (^. These energy states are not considered 
here in detail as they are irrelevant for the behaviour of the algebraically decaying 
correlations. The integral equations defining these excitations as well as analytic 
results for the gaps in various limits of p and rj can be found in the appendix. 



3 Critical exponents of the correlation functions 
and masses of charge carriers 

In this section we present the results for the correlation functions which are obtained 
from applications of conformal field theory to finite-size calculations performed as in 
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[ 17]. In the anisotropic model we have only one gapless excitation in contrast to 
the isotropic model with 25" + 1 such modes. As a result our model is described by 
a c = 1 conformal field theory and the scaling dimensions of the primary fields are 
uniquely expressed in terms of the dressed charge [ 18- 21] as 



X 



AM 



+ [avo)d]\ (12) 



where AM and d label the quantum numbers which specify the massless excitations. 
AM is the change in the number of complexes (^ compared to the groundstate, d 
is the number of complexes excited from the left Fermi point to the right one. All 
other excitations correspond to broken complexes (|^) and have a gap. Therefore 
these excitations do not affect the critical properties. 

The dressed charge function ^{v) is given by the solution of the integral equation 
associated with (H) 

p2lT — VQ 

27TC{v) = 2n- e\v -w)i{w)dw. (13) 

According to conformal field theory the long-distance asymptotics of the equal- 
time correlation functions of primary scaling fields $(r) are determined by the critical 
exponents 2x 

(<f W$(o)) ^ kF = ^/^^- (14) 

The dimensions of descendent fields differ from x by integers N^. 

We first consider the density- density correlations. The leading contributions to 
the algebraic decay of this correlator are given by formula ( P^ for d = 0, {N~^, N~) = 
(1, 0) or (0, 1), and d = 1, N"^ = 0, respectively. This leads to the asymptotic form 

(p(r)p(O)) ~ p2 + Air-^ + Aar-" cos(2A;irr), (15) 

where 

25+1 

pW= E ^rs, a = 2[^{vo)]\ (16) 

s=l 

Turning to the superconducting correlation functions we discuss the correlations 
of complexes of 2S+1 particles 

2S+1 

C{r) = n Crs. (17) 

s=l 

In this case the operator changes the number of complexes by AM = 1 with d = 
or (i = 1/2 for half- integer or integer spin S, respectively. We thus obtain 

{C+{r)C{0)) - Br-f^, p = -, for half-integer S, (18) 

a 
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and 

(C+(r)Cm - Br-'^ cosikpr), (3=- + ^, for integer S. (19) 

a 4 

In our one-dimensfonal system we have no superconductivity in the hteral sense. 
Our model does not have finite off-diagonal long-range order. However, the super- 
conducting correlations have a longer range than the density-density correlations 
provided P < a. Analytically we find a{p = 0) = 2 and a(pniax = 1) = 2/(25* -|- 1)^. 
This implies that for all S there is a density regime [0, pd where the system has dom- 
inating superconducting correlations, see Fig. 1. This property of the anisotropic 
model differs from the behaviour of the isotropic t — J model where spin and charge 
rapidities are also classified as bound states, however, with arbitrarily small binding 
energy. Therefore, the latter system does not manifest superconducting properties. 

The formation of complexes of 25+1 bound particles can be substantiated through 
a study of the conduction properties, notably the transport masses of the charge 
carriers. Following [ 22, 23,17] we determine the charge stiffness (Drude weight) 
Dc from the dressed charge ^ 

= ^^^^v^eivo). (20) 
Ztt 

vp is the velocity of the gap less excitations 

27ro-(i;o) 



and e{v) is the dressed energy defined in (|A.1|) in the appendix. Figs. 2 and 3 show 



numerical results for the charge stiffness and for the effective mass m of the charge 
carriers per bare mass rrie of the particles 

- = ^, (22) 

where D'^ = sinnj^^ is the charge stiffness of free fermions with spin S. In the 
strong-coupling limit the charge stiffness is of order e"^'^'' 

and therefore the transport masses are of order e^'^''. The reason for these large 
masses is the binding of the particles in complexes. The hopping of particles from 
one lattice site to another is always accompanied by the breaking of the complex to 
which it belongs, at the expense of the binding energy of order e^'^^. 

In summary, we have obtained the spectral gaps, the transport masses, and the ex- 
act critical exponents of the long-distance asymptotics of the correlation functions of 
an anisotropic multicomponent t — J model. In particular the algebraically decaying 
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density-density and superconducting correlations have been analyzed and supercon- 
ducting properties have been found. The one-particle Green's functions (c^^cqs) on 
the other hand show exponential decay, because the fermion operators Crs change the 
total number of particles by 1 corresponding to an excitation with gap. As a con- 
sequence the distribution function (n^^) is non-singular. Therefore the model under 
consideration is characterized as a (marginal) Luttinger liquid with a single mode of 
gapless charge excitations and 25* spin excitations with gap. 



Appendix 

Here we give a brief account of the elementary excitations over the ground state 
which is characterized by a distribution of the parameters in strings of length 
25 + 1. 

The first kind of excitations consists of a redistribution of these strings and is 
of particle-hole type. The relevant equation governing the distribution of the string 
parameters Va is (^) from which we derive the integral equation for the dressed energy 
e{v) in standard way 

^M + T^/ e'{v-w)e{w)dw = e{v) + {2S + l)fi, (A.l) 

where e{v) is the bare energy of (25 + l)-strings with rapidity v and we have intro- 
duced the chemical potential fi per particle. The chemical potential is a function of 
the particle density p (= p(fo)) via the condition 

e{vo) = e{27f - vo) = 0. (A.2) 

The dressed energy e{v) is negative for v inside the interval [vq, 27r— fo] and positive for 
V outside, thus corresponding to excitations of hole- and particle-type, respectively. 
The general particle-hole type excitations read 



E-Eo = E^(^p)-E^( 



Vj, Vn 

where the momentum 7r(f ) is defined by 

'k{v) = \— r \e{w) + {2S + l){2cos\iT]- p)i{w)]dw. (A.4) 

2 sinh rf Jn 

There are 25* massive modes of excitations corresponding to the breaking of some 
strings (^. In the general situation the Xf^ are distributed in additional strings of 
length m = 1, 2, 25. Instead of (H) we obtain 

/ 2 S* -I- 1 \ ^"^ 
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where the parameters f^"*"* are the rapidities of the strings of length m and we have 
used the definition 

(25-l+m)/2 

0m{v)= E Q{vJv)- (A.6) 

«=(2S+l-m)/2 

Applying standard transformations to (|A.5|) we arrive at the energy-momentum dis- 
persions of the massive modes in terms of the dressed energy e{v) 

, , ^ / sinh mri sinh ri , \ 1 /-STr-uo ^ , ^ -, 

em{v) = 2[ — m cosh 7] \ + mjj, / u [v — w)e{w)aw, 

\ cosh mr] — cos v j 271 Jvo 

Pm{v) = n -e(^v,^r] 

+ - — —— / Omiv - w) [e{w) + {2S + l){2coshr] - fi)^{w)]dw. 
Air smh rj Jvo 

(A.7) 

The gap = 6^(71") of the m^^ excitation is found to be and [1 — m/ (25*+ l)]e'' 
in the limits of the anisotropy parameter 77 — and 77 00, respectively. For density 
p ^ we find 

A„ = 2 sinh , ( - " ^ 1 ) ■ 

ycoshm?7 + 1 25 + 1 cosh(25 + 1)1] + 1/ 

and for maximum density p ^ 1 

sinh mr] m , _(2s+i),?fc sinhmT^fc ^ 



cosh mr] + 1 25* + 1 ^ sinh(2S' + l)r/A; y ' 

(A.9) 

in the latter limit the dressed energy is 

e{v)=4 sinh r] e-"'' . , ^'""^"^^ , (cos kv-1). (A. 10) 



fc=i 



sinh(25 + l)r]A; 
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Captions 



Figure 1: (a) Depiction of the critical exponent o; for 5* = 1/2 and different values off] — 

0. 1, 0.5, 1, 2, 10. For values of a above the broken line we have dominating 
superconducting correlations, (b) The same for 5" = 1. 

Figure 2: (a) The charge stiffness Dc ioi S = 1/2 and different anisotropics ?7 = 0.1, 0.5, 

1, 2, 5. (b) The same for S — 1. Note that the values for 77 = 5 are zero within 
the graphical resolution. 

Figure 3: (a) The effective mass m of the charge carriers for S = 1/2 and different 
anisotropies rj = 0.1, 0.5, 1, 2. (b) The same for S — 1. 
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